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a^ ; ABSTRACT 

^ : 

Ph! We investigate the role of large amplitude sub-critical thermal fluctuations 

p • in the dynamics of first order phase transitions. In particular, we obtain 

D i a kinetic equation for the number density of sub-critical fluctuations of the 

1^ • broken-symmetric phase within the symmetric phase, modeled as spherical 

^ i bubbles, and solve it analytically for temperatures above the critical tem- 

' perature. We study the approach to equilibrium and obtain the equilibrium 

rS I distribution of sub-critical bubbles of the unstable phase by examining three 

' possible mechanisms responsible for their removal; their shrinking, their cou- 

pling to thermal noise, and by thermal fluctuations of the true vacuum inside 
them. We show that for sufficiently strong transitions, either the shrinking 
or the coupling to thermal noise dominate the dynamics. As the strength of 
the transition weakens we show that sub-critical fluctuations become progres- 
sively more important, as a larger fraction of the total volume is occupied by 
the broken-symmetric phase, until the point where our analytical approach 
breaks down. Our investigation suggests that pre-transitional phenomena 
may considerably change the dynamics of sufficiently weak flrst-order transi- 
tions. We apply our results to the standard electroweak transition. 
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1. INTRODUCTION 



The study of the kinetics of first order phase transitions is by no means a new topic. 
Since the pioneer work of Becker and Doring, the condensation of supersatured vapor 
has been understood to occur by the thermal nucleation of droplets of the liquid phase 
which, when larger than a critical size, will grow and coalesce completing the transition.S 
Despite all the effort dedicated to the study of the kinetics of phase transitions, their 
intrinsic non-linear nature makes progress slow and restricted to simple systems. Even for 
simple systems in the same universality class of the Ising model many questions remain 
unanswered, from a microscopic understanding of coarse-graining to the description of 
the late stages of the transition.@ 

Our main interest is in the study of cosmological first order transitions, where the 
cooling is enforced adiabatically by the expanding Universe.! In particular, we are only 
interested in transitions where the expansion rate of the Universe is much slower than 
typical fluctuation time-scales in the system, i.e., when Mpi ^ m, where Mpi is the 
Planck mass, and m is the characteristic mass scale in the model. 

In this work we will study analytically the early stages of a first order phase transition. 
That is, we will restrict our study to temperatures above the critical temperature Tc, 
which is defined as the temperature at which the homogeneous part of the coarse-grained 
free energy density (the effective potential to some order in perturbation theory) exhibits 
two degenerate minima. By focusing on temperatures above the critical temperature, 
we will be examining the behavior of thermal fluctuations around the high temperature 
equilibrium state of the system. In particular, we will consider models described by a 
potential V{(j), T) which at high temperatures has a minimum at (0) = and that, as the 
temperature drops to some value Ti, develop a new minimum at (0) = 0_|_(T). At the 
critical temperature < Ti, V^(0, T^) = V"(0+, Tc). We will examine thermal fluctuations 

1 



around {(f)) = as the temperature drops below Ti. 

The reader may be wondering why should anyone be interested in studying the kinetics 
of first order phase transitions above Tc. After all, it is a well-known fact that first order 
transitions evolve through the nucleation and subsequent percolation of bubbles larger 
than a critical size as the temperature drops below T^. This is indeed the case for 
sufficiently strong transitions. However, for weak enough transitions the critical bubble 
picture must be modified, as has been observed for several condensed matter systems in 
the past two decades. A particularly striking example is the isotropic-nematic transition 
in certain liquid crystals; the transition is first order (there is a discontinuous jump of the 
order parameter) even though there is no release of latent heat.0 In fact, large amplitude 
fiuctuations of the nematic phase within the isotropic phase have been observed above the 
critical temperature for the transition, a typical example of the so-called pre-transitional 
phenomena in condensed matter literature.@ (Interestingly enough, the free energy density 
used to describe the isotropic-nematic transition, the Landau-de Gennes free energy, is 
formally identical to the effective potential obtained for the electroweak transition.) Here, 
we would like to investigate this possibility within the context of hot field theories. The 
point is that the usual vacuum decay mechanism for first order transitions relies on results 
obtained within homogeneous nucleation calculations, which assume that below the 
system has small fiuctuations about the metastable phase. In other words, it is assumed 
that the path integral controlling the transition rate is dominated by its saddle point, 
given by the solution to the Euclidean equations of motion. Small fluctuations are then 
incorporated by evaluating the path integral by a Gaussian approximation. But as the 
transition grows weaker, large amplitude fluctuations about equilibrium will become more 
probable and the approximations used will break down.i Instead of expanding about a 
homogeneous metastable background, one should be expanding about an inhomogeneous 
background consisting of these large amplitude sub-critical fluctuations. This is clearly 
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a very hard task which nevertheless must be undertook if we are to better understand 
the dynamics of weak first order transitions. The present work is but a first step in this 
direction as we attempt to obtain the equihbrium distribution of sub-critical bubbles 
at temperatures above Tc. This way we can obtain the fraction of the total volume 
which is occupied by the broken symmetric phase as is approached from above, and 
thus examine the validity of the homogeneous nucleation picture as a function of the 
strength of the phase transition. Gleiser, Kolb and Watkins (GKW) have first suggested 
this possibility, concentrating on correlation volume bubbles, but they did not take into 
account that subcritical bubbles are unstable and, therefore, shrink. (See, however, their 
comments in Section 3.3). Later, Anderson i remarked that subcritical bubbles may 
disappear not only by shrinking but also due to thermal noise. As a result of any of these 
two mechanisms he concluded that the fraction of the total volume filled by bubbles of 
the unstable phase will be much smaller than proposed by GKW. Here we incorporate 
these two disappearance mechanisms into a kinetic description of the transition in order 
to investigate the importance of sub-critical bubbles near the critical temperature. We 
find that sub-critical fluctuations become progressively more important as the strength of 
the transition weakens, until the point where the approximations we use in order to treat 
the problem analytically break down and our analysis cannot be trusted quantitatively. 
However, we believe that our results are indicative of the relevance of pre-transitional 
phenomena in the study of sufficiently weak phase transitions in hot field theories. 

In the next Section we will obtain the equation governing the kinetics of sub-critical 
bubbles for temperatures just below Ti. In Section 3 we obtain the expressions for the 
thermal nucleation rates needed to solve the kinetic equation. In Section 4 we solve 
the kinetic equation in three regimes; assuming that the shrinking dominates the dis- 
appearance of subcritical bubbles, assuming that thermal noise dominates, and finally 
neglecting both shrinking and thermal noise. Starting with the whole volume occupied 
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only by the stable phase, we obtain the equilibration time scale and the equilibrium 
number density of bubbles of the unstable phase in the three regimes at the particular 
temperature considered, and establish the conditions for each of the three processes to 
dominate the kinetics. In Section 5 we compare the three time scales in the context of 
the standard electroweak model. As expected, for strong enough transitions shrinking 
or thermal noise dominate and the sub-critical bubbles play a negligible role during the 
phase transition. As the transition weakens we find that a larger fraction of the total 
volume is occupied by the broken-symmetric phase. The approximations we used to 
analytically solve the kinetic equation break down when a large fraction of the volume 
is occupied by sub-critical bubbles of the broken-symmetric phase, and we cannot carry 
our study into the limit of very weak transitions. But our results clearly suggest that 
for weak enough transitions a departure of the usual vacuum decay mechanism is to be 
expected. Concluding remarks are presented in Section 6. 



2. KINETIC EQUATION 

Let n{R, t) be the number density of sub-critical bubbles of radius R at time t. Thus, 
^"^drf' number per unit volume of bubbles of radii between R and R + dR. Since 

the bubbles can shrink, their radius i? is a function of time, R{t). We will only consider 
bubbles with R> ^, where ^ is the correlation length for fluctuations around equilibrium, 
given by 

r\T) = V"{(P = 0,T). (2.1) 

Bubbles with i? ~ ^ will be statistically dominant since any larger fluctuation has larger 
free-energy and is exponentially suppressed. (Recall that for T > the free energy of 
bubbles is a monotonically increasing function of -R.i) 
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Consider a large volume V filled by the stable phase (0) = 0, which is cooled down 
from high temperatures to a temperature just below Ti . Bubbles with radius R of the new 
phase with (0) = 0+ will be thermally nucleated in the background of the phase (0) = 
with number density n{R, t); the bubbles are energetically unfavored and will shrink away 
with some velocity dR/dt. Shrinking will always be present, unless there is a stabilizing 
mechanism for the bubbles, as in non-topological solitons. i [We would like to stress 
though that not much is known about nonlinear bubble collapse. Naively, one would 
expect small, unstable bubbles to collapse in a time ~ R. In fact, recent results on the 
evolution of unstable bubbles found a remarkable enhancement of the lifetime of bubbles 
by 3 to 4 orders of magnitude, as long as their initial radius is larger than about 2.5^ 
and their amplitude at the center probes the nonlinearity of the potential.0] Thermal 
noise may also destroy small sub-critical bubbles. The importance of this effect will 
depend both on the ratio of bubble size to thermal length (i?/T~^), and on the strength 
of the coupling of the bubbles with the thermal background, which we parameterize 
with a dimensionless coefficient a. Anderson! wrote this "thermal destruction" rate as 
~ aT{RT)^. This expression would imply that the rate of disappearance of a bubble due 
to thermal noise increases with its size, a result we find counter-intuitive. Instead, we 
will conservatively take this rate to be of order aT, and thus independent of bubble size. 
As we are mostly interested in small bubbles here, we expect that this overestimate of 
the true rate will not compromise our results. In any case, it will only make our final 
results stronger, as the final fraction of volume occupied by sub-critical bubbles will be 
larger than what we obtain. This expression yields a lower rate than that of Ref. |^, since 
T > R^^. We take, therefore, the rate per unit volume of disappearance of sub-critical 
bubbles due to thermal noise to be Ttn{R) — aT/^nR^. 

We now proceed to obtain the rate equation. We can say, quite generally, that the 
number of bubbles that at time t + dt have radius R + dR in a volume ^ ^ is equal 
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to the number of bubbles at time t with radius R plus the net change in the number 
of bubbles in the range {R, R + dR) due to thermal nucleation, disappearance due to 
thermal noise and shrinking in the time interval dt, 

n{R + dR,t + dt)V = n{R, t)V + (^] To^+{R)Vdt 



V 



^ ^ T^_^oiR)Vdt - i^-^j TTNiR)Vdt. (2.2) 

Thus, to first order in dR and dt we obtain, after subtracting n{R + dR, t) to both 
terms and dividing by Vdt, 

dn{R,t) _ dn{R,t) f dR\ ^ ^^^^^^^ 



^^^T^^o{R)-(^]Ttn{R) (2.3) 



dt dR \dt J W 

Here, ro^+(-R) {T^^q{R)) is the rate per unit volume for the thermal nucleation of a 
bubble of radius R of phase = 0+ within the phase = (phase = within 
the phase 0+). These were the only two rates considered in GKW.I The volume ratios 
Vo(+)/l^ take into account the fact that the total volume in each phase changes in time 
due to the evolution of n{R,t). The initial conditions we choose are 

Vo{t = 0) = V and V+{t = 0) = , (2.4) 

that is, all volume V is initially in the phase = 0. Also, V+ must be understood as 
the volume of the (+)-phase in bubbles of radius R only, since we are following the 
evolution of n{R,t). Thus = ^iiR^n{R,t). In Eq. (|2.3| ) bubbles of radius R can 

disappear due to their shrinking (accounted for by the first term in the right-hand side 
with dR/dt < 0), due to thermal noise (last term in the right-hand side), and due to 
the nucleation of bubbles of the (O)-phase in their interior. This latter process will, in 
general, involve bubbles of (O)-phase of different radii as we will discuss later on. 
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Because 

Eq. ( p^.3| ) is an integro-differential equation and must be solved within certain approxima- 
tions, as with any Boltzmann-hke equation. In a complete treatment of the kinetics we 
should also write an equation for the fluctuating regions of the (O)-phase within the (+)- 
phase. Furthermore, there would be two contributions to Vq, a connected background 
volume and a bubble-like, disconnected volume. The topology of the two-phase system 
will, in general, be very complicated. Thus, we choose the range of temperatures in which 
it is consistent to restrict the (O)-phase to the background, as will be clear later. 

Before we move on, we comment on two other possible contributions to n{R, t) which 
we will not discuss here; a) Induced nucleation: small bubbles may act as seeds for the 
nucleation of other bubbles in their neighborhood due to the gain in surface energy. There 
should be an enhancement of the nucleation rate due to the presence of small bubbles in 
the background, very much like the presence of impurities in condensed matter systems, 
b) Collision of bubbles: small bubbles may acquire a thermal velocity and collide forming 
larger bubbles. Of course, for very short lived bubbles none of these processes should be 
very important. However, we again stress that large enough bubbles can be quite long- 
lived, and that a truly realistic scenario of weak first order transitions will be much more 
complicated than our simple model. The worth of the present effort relies on it being the 
first attempt to go beyond GKW, by incorporating more complicated out-of-equilibrium 
processes in the usual description of phase transitions in field theories. 

3. THERMAL NUCLEATION RATES 

In order to solve Eq. ( p.3| ) we need to determine the thermal nucleation rates ro(+)^+{o) • 
For temperatures T < Ti, it is reasonable to assume that most bubbles will be nucleated 

7 



4 dniR^ 



(2.5) 



with radius R > ^, where ^ was defined in Eq. (|2.1|) . For simphcity we will assume 
the same correlation length for fiuctuations around both minima, even though at Ti the 
correlations about 0+ diverge. Following GKW we write for the rates, 



0{+)-+(0) 



(R) = AT'exp -F (0+(o)) /T 



(3.1) 



where A is a constant of order unity and the ansatz for the sub-critical bubble configu- 
rations (t> is 



0+(r) = (/)+exp —r'^/R^ 



6o(r) = — exp —r^jR 



(3.2) 



with the free energy functional -F(0) of the bubble configuration given by 

r. / 7T\ 2 

F (0) = 47r y" r^dr 



1 /cf0 

2 \^ 



(3.3) 



In order to move on, we will choose a specific (but quite general) potential, 

1/(0, T) = - 7(r)03 + , (3.4) 



2 ' 4 ' ' 

where 7(T) and A(T) are positive definite functions of T and rn^iT) can be negative 
below a certain temperature T2 < T^.. This potential has a minimum at = as long as 



m2(T) > 0, and a local minimum at 0+ = [37(r) + ^j9-f^{T) -Am\T)X{T) ] /2A(T), 
which appears at a temperature Ti given by the solution of -f^{Ti) = Am^{Ti)X{Ti)/9, 
when 0+(Ti) = 37(Ti)/2A(Ti). Below this minimum becomes the global minimum 
and the minimum at = becomes metastable. Using Eqs. ( ^.2| ) and ( ^.41) in Eq. (P^), 
we obtain for the free energy F(0),0 

f(0o(+)) =ao(+)i? + /9o(+)i?' , (3.5) 

where 



3^2 



(3.6) 



/?0 



+XiT)<Pl 



-1 + 



(V2-8)+7(T)0+ (^3 

3v/2 73 1 ^ 
~8 9~ ^ 32 



3v^ VS' 



(3.7) 



and 



8 -^^(^)^^ 



32 



(3.8) 



4. SOLVING THE KINETIC EQUATION 



We will solve Eq. ( |2.3| ) in different regimes assuming that the dominant process for the 
disappearance of bubbles is: 1.) their shrinking, 2.) thermal noise and 3.) the nucleation 
of bubbles of the true vacuum inside them. This will allow us to make a direct comparison 
of the results for the relaxation time scale and equilibrium number density obtained in 
these regimes, in the context of the electroweak transition in the next Section. 

4.1 Kinetic Equation with Shrinking 

We include only the shrinking term, neglecting the last two terms of Eq. ( p.3|) , and 
solve this equation with the approximation. 



4 dn{R,t) 

TTiX 



= 1 - |/| ~ 1 . (4.1) 



3 dR 

In words, we will solve the kinetic equation in the regime in which the volume in the 
(+)-phase, given by the sum over all bubbles with radii from ,^ to oo (actually, to ~ 
\/i/3 ^ jg small, i.e. |/| <^ 1. We will check our results for consistency after obtaining 
n{R,t). With this approximation the kinetic equation becomes, 
dn{R,t) dn{R,t) 



dt dR 



[f{R)]+To^4R) , (4.2) 



where we wrote dR/dt = —f{R) < 0. The positive definite function f{R) encompasses 
the physically interesting cases dR/dt = v = constant, and dR/dt ~ We enforced 

the shrinking explicitly by choosing the minus sign in dR/dt. The rate ro-^+ was defined 
in Eq. (|3.1|) , with the free energy given by Eq. ( p.5|) with a+ and /?+ given by Eqs. (|3.6|) 
and ( |3.8| ) respectively. In order to solve this equation analytically we make one further 
approximation; we neglect the volume term in the expression for the free energy, so that 
the nucleation rate can be written as 

To-.+ {R) - AT^exp [-a+R/T] = g{R) . (4.3) 

Here we have defined the function g{R). Given that most bubbles have radii R > ^ 
this should be a good approximation that must be tested in each application. For the 
electroweak case, independently of the parameters of the model, we obtain ~ 
6.65 for T = Ti. The ratio increases for lower temperatures. 

In order to solve the above equation, first note that the time independent equilibrium 
number density at the temperature chosen, n{R), is the solution of 
dn{R) _ g{R) 



OR fiR) ' 
which is easily obtained as 



(4.4) 



n{R) = n{0 - ^^dR . (4.5) 
Choosing f{R) = v and the definition of g{R) above we obtain, 

n{R) = exp [-a+R/T] . (4.6) 



Using Eq. (|4.4| ), we can rewrite the kinetic equation as 

3Yi^ . ^/(«). (4.T) 
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where we introduced the departure from equihbrium, Y{R,t) = n{R,t) — n{R). This 
equation is solved by writing 



Y{R,t) = X(i?)exp [-t/ri] 



(4.8) 



so that Ti is the relaxation time and the solution approaches the equilibrium distribution 
as t ^ oo. Choosing again f{R) = v we obtain upon substitution in Eq. (|4.7|), 



X{R) = Xoexp [-R/vn] , 



(4.9) 



where Xq is a constant, which, together with ri, will be fixed by the initial condition 
Y{R, 0) = -n{R). Using Eq. (O) we find 



X, 







and Ti 



T 



Putting the results together we find for the solution of the kinetic equation. 



(4.10) 



n{R,t) =n{R){l-e 



-a+R/T 



1 



-a+vt/T 



(4.11) 



Thus, the relaxation time to approach the equilibrium distribution n{R) is the time for 
a bubble of characteristic radius vti = T/a+ to shrink with constant velocity v. (The 
reason why ti is independent of R is due to our choice of constant shrinking speed for 
all bubbles.) Substituting this solution into the consistency condition given in Eq. (|4.1| ) 
we obtain 



1 + 3 



T 



6 



T 



6 



T 



< 1 . 



(4.12) 



The consistency condition is then an expression of the validity of the semi-classical ap- 
proximation, which requires that exp —F (0+) /T ~ exp [—a^R/T] ^ 1; if the smallest 
bubble has a large nucleation barrier (= F^/T) its production is exponentially suppressed 
and only a negligible fraction of the total volume will be occupied by the (+)-phase. 



4.2 Kinetic Equation with Thermal Destruction 
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In Ref. 1^ Anderson writes the thermal destruction rate per unit volume as 

T{R)^^ ^ . (4.13) 

With this expression, the rate of disappearance of a bubble increases with its volume. As 
mentioned above, we think this is counter-intuitive, since larger bubbles should be less 
prone to be detroyed by thermal noise than smaller ones, not the contrary. In lack of a 
better understood expression we will take this rate to be aT, on dimensional grounds, 
independently of the size of the bubble. The rate per unit volume of the disappearance 
of bubbles due to thermal noise is then, 



m)TN = 4^ • (4.14) 



3 

The constant a is proportional to the coupling of the bubble to the thermal bath to some 
power. This is a difficult quantity to obtain microscopically, and will depend on how the 
field couples to itself and to other fields in the model. There has been some progress 
recently in the computation of the viscosity coefficient in hot field theories, although the 
results are model dependent. @ However, it is by now clear that viscosity (which is related 
to the coupling to the bath by the fluctuation-dissipation theorem) in field theories is at 
least a two- loop effect. 

The kinetic equation is now, 

dn{R,t) fVo\^ 4_ . 



^y)To-.+ {R)-^^RMR,t)^TN{R) . (4.15) 

We will again make the approximation of Eq. (4.1). Using Eqs. (4.3) and (4.13), the 
kinetic equation becomes 

^^^^ = g{R)-aTn{R,t) . (4.16) 

Thus, the time independent equilibrium number density is now 

n{R) = ^ = ^e--^/^ , (4.17) 
al a 
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and we solve the equation for (n — n) with the ansatz 

n{R, t) - n{R) = X{R)e~''^^ , (4.18) 

and the initial condition n(-R, 0) = 0, which yields X{R) = n{R). Thus, = aT and 

n(R,t) = n(R)(l - e""^*) = ll^e-"+^/^(l - e""^*) (4.19) 

a 

The relaxation time is thus T2 = {aT)~^ as expected. Comparing this value with the 
relaxation time when shrinking alone is considered ri ~ T/a+v, as obtained in the 
previous section, we see that, under the approximations we choosed, thermal noise will 
be dominant, i.e. T2 < ti, for 

a > ^ ^ (4.20) 

Here we have taken a+ ~ (f)\ from Eq. (3.6). This is an important condition, since, as we 
mentioned above, a is proportional to some power of the coupling of the field with the 
thermal bath, and 0+/T is a measure of the strength of the transition. (For successful 
baryogenesis during the electroweak scale one needs roug hly 0+/T > l.El) Thus, thermal 
noise only dominates the dynamics over shrinking for sufficiently weak transitions, and 
for strong enough coupling to the thermal bath. The consistency condition in Eq. (4.1) 
is again given by Eq. (4.12), but in this case n{C) is given by Eq. (4.17). Thus, this 
condition is again an expression of the validity of the semiclassical approximation, as 
explained below Eq. (4.17). 

4.3 Kinetic Equation Without Shrinking and Thermal Destruction 

Neglecting the shrinking of bubbles and their coupling to thermal noise, the kinetic 
equation is 

a"W')-^l^)r„_,m-f^V.-o . (4.21) 



dt W) ' W 
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where the nucleation rates have been defined before. In order to be consistent with our 

previous approach, we will again take VqIV — 1. Since {V-^/V)T^ >o{R) is the rate at 

which bubbles of (+)-phase of radius R disappear due to the nucleation of a region of (0)- 
phase in its interior, as mentioned before, V+ is the total volume of all bubbles of radius 
R (and not the total volume of (+)-phase). Thus we have V^/V = {4:7!-/3)R^n{R,t). 
Consistency with the assumption Vq/V ^ 1 as t ^ oo requires that {Att /3)R^n{R) <^ 1. 
Since most bubbles of the (+)-phase will have i? ^ i^, this condition is essentially the 
same as that of Eq. (|4.12| ), as it should. 

In principle, to obtain r_| ,o we should sum over bubbles of the (O)-phase of different 

radii which can trigger the disappearance of regions of the (+)-phase. This sum should 
be dominated by the contribution of bubbles of radius C,- In this case the dominant 
process for wiping out regions of (+)-phase is nucleation of correlation volume regions 
of the (O)-phase, as was explicitly assumed in the work of GKW. However, whenever we 

will need to write r_| >o explicitly to obtain analytic expressions, we will approximate it 

with the rate of thermal nucleation of bubbles of radius R, F+^q — ^+^o{R)- We expect 
this approximation to be good since most bubbles of (+)-phase have radii R not much 
larger than ^ anyway. The kinetic equation is, therefore, 

= To^^iR) - fRMR, t)T^^o . (4.22) 
We can rewrite it as, 

= -q{R)z{R, t), z{R, t) ^ p{R) - q{R)n{R, t) , (4.23) 
where we defined 

p{R) = To^+{R) and q{R) = —R'T+_^o ■ (4.24) 
The solution to Eq. ([4.23|) is simply z{R, t) = zoexp[—q{R)t], or, using the initial condition 
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n{R,0) = 0, 

n{R, t) = n{R) 



■ n{R) 



p{R) 



q{R) 

For a given radius it!, the relaxation time is then 

"47r 



(4.25) 



r,{R) = [q{R)]-' = 



(4.26) 



that is, the time scale for the nucleation of a bubble of the (O)-phase inside a bubble of 
radius R of the (+)-phase. Note that T3 depends exponentially on R^; bubbles of larger 
radii will approach the equilibrium distribution at a much slower rate than smaller bub- 
bles. The fastest bubbles to equilibrate are those with ~ ^. With the approximation 
r+_»o — r+-+o(-R) and using the expressions of Section 2 for the nucleation rates, the 
equilibrium distribution is given by 



47r 



R^n{R) = exp 



{P+ - Po) R' 
T 



(4.27) 



The consistency condition then states that the volume occupied by the (-l-)-phase is 
small. 



n{^) = exp 



(/9+ - (3o) e 

T 



< 1 



(4.28) 



5. APPLICATION TO THE ELECTROWEAK TRANSITION 

For a given model with a first order transition which of the approaches above better 
describes the dynamics of thermal fluctuations? We are limited to study the transition 
at T < Ti when only small sub-critical fluctuations exist. However, even for this limited 
temperature range we should be able to examine the importance of incorporating the 
shrinking of the bubbles and their coupling to thermal noise into the description of the 
kinetics. In this Section we will do this in the context of the standard electroweak model 
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using the 1-loop approximation to the effective potential. Even though recent work has 
shown that the 1-loop approximation is not adequate due to large infrared corrections 
at Tc,0 we will take the 1-loop potential as an example of how to apply our methods. 
Application to other potentials is quite straightforward. 

We will compare the different relaxation time scales obtained above, ri, T2 and r^, 
and investigate which mechanism for the suppression of regions of broken-symmetric 
phase will dominate as the strength of the transition changes. Our strategy is as follows. 
First we compare ri with r2, that is, we compare the relaxation time scale incorporating 
only shrinking (Section 4.1) with the relaxation time scale incorporating only thermal 
noise (Section 4.2). We have already shown under our assumptions that thermal noise 
dominates over shrinking if the parameter a satisfies the inequality a > a^v/T"^, Eq. 
(4.20). Since a is a free parameter in our model we have two possibilities, depending on 
the value of a. If a does not satisfy the inequality, shrinking dominates and thermal noise 
is a sub-dominant mechanism during the approach to equilibrium. Otherwise, shrinking 
is the sub-dominant mechanism. Since we then will know the values of a for which 
shrinking or thermal noise dominates as a function of the strength of the transition, 
we can proceed by comparing these two processes with the inverse thermal nucleation 
rate time scale, given by (Section 4.3). As a result, we will be able to establish which 
mechanism dominates the dynamics as the transition's strength varies for different values 
of a. 

We take the Higgs mass to be the parameter that controls the strength of the transi- 
tion, while we fix the top mass at mx = 130 GeV. For the 1-loop approximation to the 
electroweak potential, we can writeEl 



where the constants D and E are given by = [6{mw/o')^ + '^{mz/cry + 6(mT/o")^] /24, 




(5.1) 
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and E = [6{mw/o')'^ + 3(m^/o")'^] /127r ~ 10 ^. T2 is the temperature at which the origin 
becomes unstable, given by 



T2 = ,J{ml-8Ba^)/AD , (5.2) 

where the physical Higgs mass is given in terms of the 1-loop corrected A as m'jj = 
{2\ + 12B)a^, with B = (6m^ + 3m| - 12m^) /647r V. We use mw = 80.6 GeV, 
rriz = 91.2 GeV, and cr = 246 GeV. The temperature corrected Higgs self-coupling is 
1 



Xt — A — ■ 



,(5.3) 



167r2 

where the sum is performed over bosons and fermions (in our case only the top quark) 
with their respective degrees of freedom gB(F)- Also, In cb = 5.41 and In cp = 2.64. 



This potential is equivalent to the potential of Eq. (3^), with the replacements 
m2(T) = 2L)(T2 - T|), 7(T) = ET, and A(T) = At. The temperatures Ti and 
are given by 



Ti = Ts/^l - gEySXrD, and = Tg/^l - E^XtD . (5.4) 

From these values for Ti and Tc, it is clear that the quantity x = E'^/XtD can be 
used as a measure of the strength of the transition. For x = 0, the transition is second 



order. As shown in Ref. [T^, x ~ 0.027 at mu = 60 GeV, becoming smaller for larger 
Higgs masses. The transition is already quite weak at the experimental lower bound of 
In Fig. 1 we show the quantity function of the Higgs mass for 

niH > 35 GeV. 

With the potential above, we can easily compute the nucleation rates using the ex- 
pressions of Section 2, as done by Gleiser and KolbJll We have then all that is needed 
to compute the relaxation time scales ri, T2, and T3(^) [since T3 depends on R, we con- 
sider only the fastest bubbles to equilibrate, those of radius ~ and the equilibrium 
distributions n{R). 
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Let us start by comparing the relaxation time scales for the shrinking term and the 
thermal destruction term, following Eq. (4.20). In Fig. 2 we plot a;+(Ti)/Tf (we took 
w = 1) as a function of the Higgs mass. Thermal noise dominates over shrinking if a lies 
above the curve, for a given Higgs mass. Notice that for small Higgs masses, i.e. strong 
transitions, shrinking should be the dominant process unless the coupling to the bath is 
unrealistically large (or non-perturbative). However, since there are uncertainties in the 
precise expression for the thermal destruction rate, we will consider both possibilities, 
with shrinking and thermal noise dominating the dynamics, respectively. If shrinking 
dominates, we can neglect thermal noise and compare ti and T3(^) directly. This is done 
in the continuous lines of Fig. 3 where we show the ratio of Ti/r3(^) as a function of 
the Higgs mass. The results are shown for two different shrinking velocities, v = 1 and 
V = 0.4. For both velocities, with Higgs masses below 52 GeV we find that ti ^ ^3(0) 
so that the approach to equilibrium is dominated by the shrinking. In Fig. 4 we show 
the two consistency conditions Eqs. ( [4.12| ) and ( [4. 281) as a function of the Higgs mass. 
We see that the fraction of the total volume occupied by the (+)-phase increases rapidly 
as the Higgs mass increases. That is, as the transition grows weaker a larger fraction 
of the volume is occupied by the broken-symmetric phase. Even though our results 
indicate that we cannot trust our approximations for uih ^ 53 GeV, where V+/V > 10%, 
it is quite clear that the reason our approximations break down is precisely the large 
fraction of the volume occupied by sub-critical bubbles. These results suggest that pre- 
transitional phenomena will be present for weak enough transitions and that we should 
expect modifications of the vacuum decay picture in this case. A precise description of 
the transition is beyond the scope of the present formalism. 

We now consider the case in which thermal noise dominates over shrinking. By 
inspecting Fig. 2, we can find an appropriate value of a which is large enough for the 
inequality Eq. (4.20) to be always satisfied. We will take a = 0.5 as an illustration, even 
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though this value is probably unrealistically large. In Fig. 3 the dashed line shows the 
ratio T2/t3(^) as a function of the Higgs mass. Note that again thermal noise becomes 
sub-dominant as the Higgs mass is increased, although it does so at a slower rate than the 
shrinking term. However, the qualitative conclusion is the same in both cases. For large 
enough Higgs masses, the fraction of the total volume in the broken-symmetric phase 
becomes substantial due to the increasing weakness of the transition. This weakness 
is characterized by the fact that both the shrinking of sub-critical bubbles as well as 
their destruction due to thermal noise become less important in the description of the 
transition. Since the breakdown of our approximations is due to the failure of the semi- 
classical, or dilute gas approach, we are confident that as the Higgs mass increases beyond 
the limit of validity of our approximations a regime will be reached in which a departure 
from the usual false vacuum decay mechanism is to be expected, although at this point 
we cannot afford to make a quantitative prediction. A rough estimate may be obtained 
by comparing, at the nucleation temperature, the typical distance between sub-critical 
bubbles to the radius of a critical bubble obtained by the usual calculations. If the 
distance between sub-critical bubbles is of the order of the critical radius, the usual 
nucleation mechanism must be revised. Work on this topic is currently in progress. 

6. CONCLUDING REMARKS 

We have obtained a kinetic equation describing the approach to equilibrium in first- 
order transitions for temperatures above the critical temperature. For sufficiently strong 
transitions and making approximations that we consider reasonable, we were able to 
analytically solve the kinetic equation for three different regimes, determined by the 
dominant mechanism responsible for the suppression of regions of the broken-symmetric 
phase within the symmetric phase. The three processes are the shrinking of the sub- 
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critical bubbles, their destruction due to the thermal bath to which they couple, and by 
thermal nucleation of regions of the symmetric phase in their interior. By obtaining the 
relaxation time-scales in all regimes, we were able to study the relative importance of 
each process in the early stages of the transition. 

We applied our approach to the standard electroweak transition, showing that for 
Higgs masses below 55 GeV or so, the approach to equilibrium is dominated by either 
shrinking or thermal destruction, depending on the strength of the coupling of bubbles to 
the thermal bath. In this case, the total volume occupied by the equilibrium distribution 
of sub-critical bubbles of the broken-symmetric phase is negligible, and the transition 
proceeds by the usual vacuum decay mechanism. As the Higgs mass increases and the 
transition becomes progressively weaker, we found that a larger fraction of the total 
volume becomes occupied by the broken-symmetric phase, forcing the breakdown of 
our analytical approximations. However, it is clear from our results that a regime will 
eventually be reached in which a substantial fraction of the volume is in the broken- 
symmetric phase as the critical temperature is reached. In this case we should expect a 
departure from the usual vacuum decay mechanism. We cannot be quantitative about 
the value of the Higgs mass (for a given top mass) where this occurs due to the limitations 
of our analytical approach. Given that the lower bound on the Higgs mass is now above 
60 GeV, we expect very interesting physics to be lurking behind our present knowledge 
of the dynamics of weak first-order transitions. 
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FIG. 1. The parameter x = E'^ / \tD as a function of the Higgs mass for several values 
of the top mass. 



FIG. 2 The ration of relaxation time-scales for shrinking and thermal destruction 
processes as a function of the Higgs mass. The region above the curve defines the domain 
for thermal destruction domination. 



FIG. 3. The continuous lines denote the ratio between the shrinking and reverse 
nucleation time-scales, as a function of the Higgs mass for shrinking velocities v = 1 
and V = 0.4. The dashed line denotes the ratio between thermal destruction and reverse 
nucleation time-scales with parameter a = 0.5. 



FIG. 4. The consistency conditions of Eqs. |4.12| and [4.28| as a function of the Higgs 
mass. 
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